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ABSTRACT 

We investigate the emission of rotating magnetized neutron stars due to the acceleration and radiation 
of particles in the relativistic wind and in the magnetotail of the star. We consider that the charged 
particles arc accelerated by driven collisionless reconnection. Outside of the light cylinder, the star's 
rotation acts to wind up the magnetic field to form a predominantly azimuthal, slowly decreasing 
with distance, magnetic field of opposite polarity on cither side of the equatorial plane normal to 
the star's rotation axis. The magnetic field annihilates across the equatorial plane with the magnetic 
energy going to accelerate the charged particles to relativistic energies. For a typical supersonically 
moving pulsar, the star's wind extends outward to the standoff distance with the interstellar medium. 
At larger distances, the power output of pulsar's wind E w of electromagnetic field and relativistic 
particles is redirected and collimated into the magnetotail of the star. In the magnetotail it is proposed 
that equipartition is reached between the magnetic energy and the relativistic particle energy. For 
such conditions, synchrotron radiation from the magnetotails may be a significant fraction of E w for 
high velocity pulsars. An equation is derived for the radius of the magnetotail r m (z') as a function 
of distance z' from the star. For large distances z' , of the order of the distance travelled by the star, 
we argue that the magnetotail has a 'trumpet' shape owing to the slowing down of the magnetotail 
flow. We compare results with the Guitar and Mouse nebula, and conclude that the tail of the Mouse 
may be connected with the long magnetotail behind the pulsar. We argue that the shock waves and 
elongated structures may also be observed in a misdirected or shutoff pulsars and may be used as a 
tool for finding these objects. 

Subject headings: stars: neutron — pulsars: general — — stars: magnetic fields — X-rays: stars 



1. INTRODUCTION 

The relativistic winds of supersonically moving pul- 
sars are observed to form shock waves in the interstellar 
medium (ISM) and in some cases prominent wakes or 
"tails" . Several pulsars are known to power bow shock 
waves and tails, the "Guitar" nebula (Cordes et al. 1993; 
Chatterjee & Cordes 2002), millisecond pulsars (Kulka- 
rni & Hester 1988; Bell et al. 1995), the "Duck" neb- 
ula (e.g., Thorsett, Brisken, & Goss 2002), the "Mouse" 
nebula, observed in the X-ray and radio bands and pow- 
ered by the spin-down of a 98 ms pulsar (Camilo et al. 
2002; Gaensler et al. 2004), the pulsar B0740-28 (Jones, 
Stappers, & Gaensler 2002), and the neutron star candi- 
date RX J1956. 5-3754 (van Kerkwijk & Kulkarni 2001a, 
b). Axisymmetric simulations of the interaction of mag- 
netized stars with the interstellar medium (ISM) have 
shown that such interaction leads to formation of long 
magnetotails behind the stars (Romanova et al. 2001; 
Toropina et al. 2001). Such magnetotails may allow rel- 
ativistic particles from the pulsar to propagate to large 
distance behind the pulsar. 

In isolated pulsars a relativistic magnetohydrodynamic 
(MHD) wind forms outside of the light cylinder and 
propagates outward (Goldreich & Julian 1969, hereafter 
GJ69; Arons & Tavani 1994; Arons 2004). This wind 
forms a magnetically dominated, azimuthally wrapped 
disk-like equatorial structure which expands with veloci- 
ties of the order of the speed of light and is driven mostly 
by the magnetic pressure. Observations reveal disk-like 
structures around a number of pulsars, including one in 
the Crab nebula (Weisskopf, et al. 2000), the Vela pul- 



sar (Pavlov et al. 2001), and a few other cases (Helfand, 
Gotthelf, k, Halpern 2001; Gotthelf 2001). The orienta- 
tion and thickness of these disk-like structures (or, PWN 
tori) were estimated by Ng and Romani (2004) who sug- 
gested that particles are accelerated more efficiently in 
these disk-like structures. It was also noticed that the 
rotation axis of the disk approximately coincides with 
the direction of pulsar's velocity, so that the geometry of 
the flow is approximately axi-symmetric (Spruit & Phin- 
ney 1998; Lai, Chernoff, & Cordes 2001). In a few cases 
jets were observed which are approximately aligned with 
the direction of motion (e.g., Weisskopf, et al. 2000). 

Formation of such an equatorial disk structure was re- 
cently modelled in axisymmetric relativistic MHD sim- 
ulations by Del Zanna, Amato, and Bucciantini (2004), 
and Komissarov and Lyubarsky (2004) for the case of an 
intial split monopole magnetic field. They have shown 
that a magnetically-dominated disk is formed around a 
rotating pulsar and that this disk has a small thickness. 
Formation of rotating, azimuthally wrapped disk-like 
structures were also observed in non-relativistic axisym- 
metric simulations of an accreting, rotating star with an 
aligned dipole magnetic field in the propeller regime (Ro- 
manova et al. 2003). This disk- like magnetically domi- 
nated bulk flow carries most of energy out of the pulsar, 
however, it may be invisible unless the particles are ac- 
celerated by some mechanism. 

Different models have been proposed to explain accel- 
eration of particles in this wind. In one class of models 
ideal relativistic magnetohydrodynamics is assumed to 
hold outside of the light cylinder with the acceleration 
of plasma due to the spatial variation of the electric and 
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magnetic fields (e.g., Vlahakis 2004). In these models 
it is suggested that the Poynting-flux to mattcr-encrgy- 
flux ratio (a) evolves from a value much larger than unity 
near the light cylinder to a value less than unity at large 
distances where the pulsar wind encounters the interstel- 
lar medium. In another class of models non-ideal plasma 
effects are considered to be responsible for the particle 
acceleration. The particle acceleration may be stochas- 
tic as discussed in the shock wave acceleration models 
(Arons & Tavani 1994; Arons 2004; Spitkovsky & Arons 
2004). Alternatively in collisionless reconnection and an- 
nihilation the magnetic field may cause random (Coroniti 
1990) or bulk acceleration of the particles (Lyubarsky & 
Kirk 2001; Kirk & Lyubarsky 2001). These models argue 
that the MHD wind consists of an azimuthally wrapped 
dipole magnetic field of the star with regions where the 
magnetic field reverses direction. These regions are likely 
sites for reconnection and particle acceleration. 

In this work we investigate a model where particles 
are accelerated in the pulsar's equatorial neutral layer 
across which the predominantly azimuthal magnetic field 
reverses direction. There may also be oppositely directed 
Poynting flux flows aligned with the star's rotation axis 
analogous to those generated from the time-dependent 
expansion of magnetic field loops of an accretion disk into 
a low density coronal plasma discussed by Lovelace and 
Romanova (2003) and simulated by Lovelace, Gandhi, 
and Romanova (2004) with a relativistic electromagnetic 
particle-in-cell method. However, note that numerical 
solutions of the "pulsar equation" for an aligned station- 
ary magnetosphere by Contopoulos, Kazanas, and Fendt 
(1999) and Gruzinov (2005) do not show Poynting flows 
along the rotation axis. 

Our model is qualitative in that we do not have global 
self-consistent calculations of the electromagnetic fields 
and sources. Compared to Coroniti (1990) we suppose 
that the misalignment angle between the rotational and 
magnetic axes is not large. Then, instead of multiple 
neutral layers in a system of magnetic stripes, a global 
extended neutral layer forms in the equatorial plane. The 
equatorial layer is not expected to be thin everywhere. 
However, the systematic acceleration will be in the radial 
direction. Thus the bulk motion of plasma and magnetic 
field will be in the radial direction mainly in the equa- 
torial plane. Thus, we expect to have a high-er plasma 
near the light cylinder which gradually changes to a much 
lower a far from the pulsar. Accelerated particles will 
gradually diffuse vertically away from the thin neutral 
layer, and their interaction with the magnetic field will 
produce the synchrotron radiation. In this paper we de- 
velop this model in semi-quantitative detail. 

The acceleration of particles in collisionless reconnec- 
tion has been discussed by a number of authors (Alfven 
1968; Dessler 1969, 1971; Speiser 1970; Cowley 1971, 
1973; Bulanov & Sasorov 1976; and Vasyliunas 1980; 
Burkhart, Drake, & Chen 1991). The present work uses 
the model by Alfven (1968) and its relativistic counter- 
part (Romanova & Lovelace 1992; Zenitani & Hoshino 
2001; Larrabee, Lovelace, & Romanova 2003). Both the 
simulations (Zenitani & Hoshino 2001) and the analytic 
theory (Larrabee et al. 2003) indicate a power law dis- 
tribution of accelerated particles. 

Most pulsars propagate through the ISM highly super- 
sonically so that the pulsar wind forms a strong shock 



wave in the ISM. A model of this interaction was first 
discussed by Schvartsman (1970). Subsequently, the 
shape of the shock wave was calculated for the hydro- 
dynamic flow by Baranov, Krasnobaev & Kulikovsky 
(1971), Baranov & Malama (1993) and by Wilkin (1996). 
However, for the pulsar case it is necessary to include the 
relativistic magnetized wind of the star. Recent models 
of the Solar wind interaction with the ISM include the 
influence of the magnetic field for both the interstellar 
and interplanetary components, and these models show 
that the shape of the bow shock may depend on these 
fields (e.g., Linde et al. 1998; Zank 1999; Bucciantini et 
al. 2004; Ophcr et al. 2004). However, only relatively 
low Mach numbers are considered and the Alfven speeds 
are typically small compared with the flow speeds. In 
the case of pulsars the Mach numbers are much larger 
than unity and the Alfven speeds may be comparable 
with the flow speeds. Strongly magnetized magnetotails 
were observed in axisymmetric simulations of rotating 
and non-rotating strongly magnetized stars (Toropina et 
al. 2001). 

In this paper we discuss the nature of the magneto- 
tails which form behind the pulsars. We argue that most 
of the power of the star's wind of relativistic particles 
and magnetic field is redirected by the internal shock 
wave and subsequently propagates down the magneto- 
tail. In several cases observations show cometary struc- 
tures aligned with the direction of propagation of the 
pulsar (e.g., Stappers et al. 2003; Gaensler et al. 2003) 
which may be the result of the pulsar wind propagating 
down the magnetotail. 

Both shocks and magnetotails may appear in the case 
of misdirected pulsars or longer-period, shutoff pulsars. 
We discuss the possible nature of the isolated neutron 
star candidate RX J1956. 5-3754 which has a shock wave. 

In §2 we discuss basic considerations. In §3 we discuss 
particle acceleration due to reconnection in the equato- 
rial neutral layer of the star, and in §4 the associated 
synchrotron radiation. In §5 we discuss the bow shock 
wave and the related magnetotail and its synchrotron 
emission. In §6 we discuss examples of pulsars with tails. 
In §7 we consider the possibility of observing wind and 
magnetotail emitting stars in the Solar neighborhood. In 
§8 we give the conclusions of this work. 

2. BASIC CONSIDERATIONS 
A pulsar's light cylinder radius is 

r t = - = ^«4.8xl0 9 Pcm, (1) 

where P is the period in seconds. The pulsar spins down 
owing to the relativistic wind which flows outward be- 
yond the light cylinder distance (GJ69). This wind con- 
sists of an electromagnetic field and relativistic particles 
and has a power output 

E w « Blr\c = ^ « 5.8 x 10 31 ^| ^ , (2) 

(GJ69), where B L = njr\ is the magnetic field strength 
at the light cylinder radius, and fj, 30 — ^/(10 30 Gcm 3 ) is 
the star's magnetic moment. The star's surface magnetic 
field is B = £?i 2 10 12 G, with B 12 = M3o/^6, where i? 6 = 
-R*/10 6 cm and the star's radius. 
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Fig. 1. — Sketch of the magnetic field configuration of misaligned rotating neutron star outside of the light cylinder, r > tl = c/ui*. For 
the case shown fl • (J, > the magnetic field is wrapped in clock-wise spiral in the upper half-space and a counter-clockwise spiral in the 
lower half-space. A neutral layer is formed in the equatorial plane where there is forced reconnection or annihilation of the magnetic field. 
The magnetic field energy goes into accelerating particles to relativistic energies. The helical field lines around the ±z axes correspond to 
possible Poynting outflows. 



Assuming the pulsar wind approximately isotropic, the 
pressure of the wind varies as p w — E w j (4irR 2 c) for dis- 
tances R > Tl from the star. For a pulsar moving su- 
personically through the ISM, the pressure of the pulsar 
wind p w (z) balances the ram pressure of the ISM pi sm v 2 
at the stagnation point at a distance z s h in front of the 
pulsar. This is the location of the bow shock of the pul- 
sar. That is, p w (z s h) = PismV 2 so that 



Z s h 



(AirpisinV 2 ) 1 / 2 



10 1 



1/2 



cm 



v 8 P 2 

(Schvartsman 1970; Baranov et al. 1971; Wilkin 1996), 
where Pism/mp ~ ni sm (in units of cm~ 3 ) is the number 
density of the ISM with m p the proton mass. The lateral 
width of the bow shock is « z s h- Thus the power which 
goes into thermal heating of the shocked interstellar gas 
is 



Esh » -zP*smV 3 TTZ 2 h » 0.9 X 10 



28 v sB 2 2 erg 



T- < 4 > 



^gftarii.- ---sri -~ p4 

where the 3/16 factor assumes a specific heat ratio of 
5/3. The power E s h is extracted from the kinetic energy 
of the motion of the star through the ISM. The slowing 
of the star's motion is however negligible. Notice that 
the ratio 



E,„ 



3 v 
64c 



1.56 x 10" 4 v 8 



(5) 



is much smaller than unity. 

3. MAGNETIC RECONNECTION 

Outside of the light cylinder the pulsar's magnetic field 
is predominantly azimuthal. If the magnetic moment of 
the star is perpendicular to its rotation axis, then the 
azimuthal magnetic field forms "stripes" of opposite po- 
larity which move outward relativistically (e.g., Michel 
1971, Coroniti 1990, Lyubarsky & Kirk 2001). Parti- 
cles may be gradually accelerated by reconnection of the 
oppositely directed fields (Coroniti 1990). Most of the 
power flow is thought to be in the form of a magneti- 
cally dominated wind. That is, the ratio of the magnetic 
energy density to the kinetic energy density of the parti- 
cles, the <j parameter, is much larger than unity. As the 
distance from the light cylinder increases, a decreases. 

We consider the case where the angle between the mag- 
netic and rotation axes is not very large. The nature of 
the envisioned driving mechanism is sketched in Figures 
1 and 2. This work focuses on the equatorial region of 
the flow, but for completeness the figures show possible 
Poynting outflows along the ±z axes. The plasma be- 
havior is non-ideal; it is not described by ideal MHD. 
We use both cylindrical (r, <p, z) and spherical (R, 9, <p) 
non-rotating coordinates to describe the electromagnetic 
field and particle motion for R > tl- The rotation of the 
star wraps the field around the rotation axis so that it is 
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Fig. 2. — The figure shows the current flows (J), electric fields (E), and Poynting vectors (S = cE X B/4-7r) of the envisioned configuration 
outside of the light cylinder, r > tl = c/uj*. 



predominantly azimuthal as shown in Figure 1. For the 
case assumed here, a;* • fi > 0, the field lines are wound 
up clockwise for z > and counter-clockwise for z < 0. 
Consequently, there is a neutral layer in equatorial plane 
where the azimuthal magnetic field changes from > 
for z < to Be/, < for z > 0. Within the neutral layer, 
the magnetic field is approximated as 

(6) 

where Az(r) is the half-thickness of the neutral layer 



which is assumed less than r, and where Bl = n/r\ w 
9.2 fiso/P 3 G. Thus the magnetic field has form of 
Archimedes' spiral with field lines given by r = const — 
tl4>. Equation (6) neglects the effect of the annihilation 
of the field. 

In the neutral layer Ampere's law gives dB^/dz = 
—4irJ r /c. Thus the change in the azimuthal field across 
the layer is 



AS, 



B d 



-Az 



2I r 
c r 



(7) 
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Fig. 3. — Sketch of the external shock, the magnetopause, and the magnctotail of a pulsar moving supcrsonically through the interstellar 
medium. The velocity of the star is assumed parallel to the star's rotation axis The solid circles indicate the shocked interstellar gas 
with sample streamlines labelled by v. The open circles indicate relativistic leptons (£) which flow into the magnetotail. The radius of the 
magnctotail is r m (z) and the bulk velocity of the magnctotail plasma is v m (z). The internal shock of the pulsar wind is indicated by "int. 
shock." The standoff distance of the shock, z s ^, is given by equation (3). 



Thus, the radial current carried by the neutral layer is where /130 = /i/(10 30 G cm 3 ) is the magnetic moment 

_ . , ..ii U30 . /„x of the star. This radial equatorial current flow acts to 

I r = B L r L c « 4.4 x 10 11 ^ A , (8) 
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magnetically pinch the particles within the neutral layer. 
An analogous current layer exists in the solar wind (e.g., 
Bertin & Coppi 1985). 

The equatorial outflow of current is balanced by an 
axial inflow of current as indicated in Figure 2. This 
current corresponds to a particle outflow rate of 



?^« 2.7 x 10 30 ^ 
dt e P 2 



(9) 



There may be oppositely directed Poynting flux out- 
flows along the rotation axis with S z = [c/ Ai:)E r B^ < 

for z < flow outward along the star's rotation axis. 
These outflows would be analogous to the Poynting out- 
flows from accretion disks (Lovelace & Romanova 2003; 
Lovelace et al. 2005). 

For the considered case f2 • fi > 0, the electric charge of 
the neutral layer is dominantly positive. Thus the charge 
of the layer gives rise to an outward radial electric field 
E r > and an axial field E z > for z > and E z < for 
z < as shown in Figure 2. The electric field gives rise 
to an axial ExB = E r B^i (< for z > and > for 
z < 0) drift of the azimuthal field and associated plasma 
into the neutral layer. This constitutes the driving force 
for the reconnection. 

Additionally, there is a radially outward E x B = 
—E z B^i drift of the particles in the neutral layer. The 
E x B drifts are of course in the same direction as the 
Poynting vector S = (c/47r)E x B which is shown in Fig- 
ure 2. 

Particles are accelerated in the electric field E r which 
decreases with the distance as 1/r or faster. This depen- 
dence is determined by the fact that the energy comes 
from the magnetic field, ~ which varies as 1/r (see 
equation 6). There is an additional factor in equation (6), 
z/Az, which takes into account the fact that at larger 
distances the relative thickness of the neutral layer will 
likely increase due to instabilities and turbulence. The 
instabilities may be driven by the misalignment of the 
rotation fl and magnetic /x axes. This is expected to 
lead to faster than 1/r decrease of the field with the dis- 
tance. In addition, part of the energy of the field will go 
to acceleration of particles, which will cause E r to fall off 
more rapidly than 1/r. Thus, we introduce a parameter- 
ized dependence of the electric field within the neutral 
layer in the form 



E r = a E B L 



(10) 



where ctE is a dimensionless quantity assumed to be a 
constant less than unity, and q is another constant which 
characterizes how fast this electric field decreases with 
the distance. The overall electrical neutrality of the equa- 
torial and axial regions implies that q > 1 . At the bound- 
ary of the neutral layer E z (r, Az) w o^Sl^l /r), where 
a' E is further dimensionless quantity assumed to be ap- 
preciably less than unity. 

Neglecting for the moment radiative energy losses, pos- 
itively charged particles drifting into the neutral layer are 
accelerated in the direction f + (r^/r)^ which is approx- 
imately radial for r > r^. The radial acceleration gives 



where m is the particle rest mass and e is its charge 
(Alfven 1968). The approximation involves assuming 
that the particles move outward with speed ~ c. In- 
tegrating eq. (11) from r = Tl to r >> r^, we find 
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eB L r L a E 



mc 2 (q — 1) 
for r S> r E and q not close to unity. 



7 w 2.6 x 10 7 



q-1 



m 



(12) 



(13) 



d_ 

TV 



(mc 2 7) w c—(mc 2 -f) w qa E cB L ( — , (11) 
ar V r / 



for the case of leptons. For protons and heavier ions the 
numerical factor is < f.4 x f0 4 . A smaller number of 
negatively charged particles is accelerated back towards 
the pulsar. 

The outward kinetic energy flux of the particles accel- 
erated in the neutral layer is 

t, mc 2 ( 7 -l)f = (^ C )^- 
at q — 1 

.MxlO-^&S (14) 

for r 3> rL, neglecting the synchrotron losses. This ki- 
netic energy flux is a fraction ctE/{q— 1) of the Goldrcich- 
Julian power (B\r\c). 

We expect different behaviors for the case where q is 
appreciably larger than unity (q — 1 ~ 1), and the case 
when q is close to unity (q — 1 <C 1). In the first case 
particles are accelerate mainly in the vicinity of the light 
cylinder. In the second case, they are accelerated along 
the neutral layer out to much larger distances. The total 
energy, however, can not be larger than the rotational 
energy of the pulsar, that is oe/i — 1 < 1 • 

Particles, accelerated in the neutral layer will be gradu- 
ally scattered by irregularities (Alfven waves) in the mag- 
netic field and will drift away from the neutral layer. At 
the same time there may be magnetic confinement of par- 
ticle to the neutral layer out to large distances. The force 
balance (or magnetic pinch) condition for the neutral 
layer follows from the momentum equation dT zz /dz = 
or T zz = const (Bertin & Coppi 1985), where Tjk is the 
stress tensor for the particles and fields. This condition 
implies that as — (Az/r)[l — (a' E ) 2 ]. The compressive 
pinch force on the neutral layer is responsible for the 
reconnection at z = being driven. 

With increasing distance from the light cylinder the ra- 
tio of the Poynting to particle energy fluxes a decreases. 
The variation of a with distance may be anisotropic, with 
smaller a in the disk and larger a above and below the 
disk. With increasing r, a larger fraction of the mag- 
netic energy is converted to particle energy in the neu- 
tral layer. Particles accelerated in the neutral layer may 
subsequently scatter so as to give a more spherical wind. 
The scattering and diffusion processes of particles accel- 
erated in the neutral layer remains to be investigated in 
detail. For simplicity of the following analysis we suppose 
that a < 1 before the bow shock is reached so that the 
flow is super-Alfvenic and supersonic at the bow shock. 

4. SYNCHROTRON RADIATION 

The acceleration of leptons in the neutral layer is in- 
hibited due to the energy lost to synchrotron radiation. 
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Including the synchrotron losses gives 



df 



To 

fq 



k"f 2 



(15) 



Here, r = r/r L , 70 = ea E B L r L /(mc 2 ), k = 
(2/3)PiiB-[rLrl/(mc 2 ), with r e = e 2 /(mc 2 ) the classi- 
cal electron radius, and [3r < 1 accounts for the fact 
that the magnetic field close to the neutral layer is less 
than Bl/t. We assume that 7 ~ 1 at the light cylin- 
der distance f = 1. The synchrotron radiation of ions is 
negligible. Equation (15) represents a simplified model 
of the synchrotron losses from the neutral layer in that it 
neglects the fact that the leptons tend to be focused to 
the z = plane where the magnetic field is appreciably 
reduced. However, this focusing will be offset if as ex- 
pected there are significant irregularities (turbulence) in 
the plasma and the magnetic field as in the Solar wind. 

For simplicity we examine the case where q = 2. It is 
clear that the Lorentz factors will not exceed 



7T7 



-( 



1/2 



/ iea E r L 



1/2 



~ 3 - lxl0 '(lf (16 > 

Thus the synchrotron losses are significant for pulsar pe- 
riods P shorter than the period which gives jmax = 7o- 
This critical period is 

(17) 

For P < P cr the synchrotron energy loss rate is 



E s 



drI r E r = ctEB 2 L r 2 L c 



(18) 



which is a fraction a^of the Goldreich- Julian power. For 
P < P cr the kinetic energy flux of the particles is 

TP -I - 11 2 dN 

-t^kin — \ fmax ±)TYIC 



7 X 10 



dt 

1/2 1/2 

31 a E \ M 3 er § 



Pr 



pl/2 



(19) 



12.2(^/^/2) < x 



For P < P cr we have E kin /E ayn 
for a E = 0.1 = (3 R . 

For P > P cr the synchrotron losses are small compared 
with power input to the particles in the neutral layer. 
The total synchrotron radiation from the neutral layer 
(r > r L ) is 



E. 



syn 



2 2 o R 2 1 dN 



«3.9xl0^(^) 2 ^^ (20) 

The kinetic energy input to the particles in the neutral 
layer is a E B 2 L r\c. 

The synchrotron radiation of the particles accelerated 
in the neutral layer is beamed radially outward in the 



equatorial plane of the pulsar. That is, the radiation is in 
a fan beam perpendicular to the pulsar's rotation axis. In 
reality, the neutral layer is probably not a thin layer with 
small width Az everywhere. The accelerated particles 
will scatter from irregularities in the neutral layer and 
will interact with the magnetic field at larger distances, 
z >> Az. This process may explain the radiation from 
the disk-like pulsar wind nebulae such as those in the 
Crab and Vela nebulae. 

For P < P cr the frequency of the peak of the radiation 
spectrum is 

1 eB Ll 2 ^tl^ 



47r mc 



1.3 xlO 22 mf^) Hz 



0R 



(21) 



where the lepton Lorentz factor is from equation (16). 
This frequency corresponds to a photon energy of w 54 
MeV. For P > P cr , the Lorentz factor of the leptons is 
70 so that 

, sy „«8.5xl0 19 (^) 2 ^Hz. (22) 

This frequency corresponds to w 350 keV. These num- 
bers are of course estimates because the value of a E is 
not known. 

5. THE MAGNETOTAIL 

Beyond the light cylinder the pulsar relativistic wind 
carries the power E w = B 2 L r 2 L c in the electromagnetic 
field and relativistic particles (GJ69). The spin down of 
the pulsar is given by Iuj^lu^ = —E w which implies that 



p2 



1+* 



_ ±± . 



8vrV 



(23) 



where t is the age of the pulsar, P is the initial period 
of the pulsar, and r w 1.4 x 10 n (P /0.02s) 2 /Aiio s ~ 
4, 300(P /0.02s) 2 /m1o Y r assuming I = 10 45 g cm 2 . After 
a time t the total energy put out by the pulsar is E w = 
Iu 2 o[t/(t + t)]/2, where c^ = 2tt/P . 

The pulsar wind also carries toroidal magnetic flux out- 
ward. The rate of transport of positive toroidal flux is 



$ + = nB L r L c « 4.1 x 10 



21 ^30 G cm 2 



P 2 



(24) 



After a time t the total (positive) toroidal flux put 
out by the pulsar is $ + = {n [iu) 2 Q T / c) ln(l + t/r) = 
(7r/2)/x/c 2 ln(l + t/r). Of course the net flux is zero. 

For early times after the pulsar's birth there is a spher- 
ical Sedov- Taylor "bubble" of plasma of radius 1Z(t) oc t a 
with a < 1. Here we consider late times in the sense that 
the pulsar has had time to move outside of this bubble. 
That is, we consider t > TZ(t)/v, where v is the pulsar's 
velocity. This corresponds to t > 3 x 10 9 s(lZ pc /v$), where 
TZp C is the bubble radius in pc, and vg = t>/10 8 cm/s. 

Figure 3 shows a sketch of the late stage of the mag- 
netotail of a neutron star moving supersonically with ve- 
locity v through the interstellar medium. For simplicity 
we have assumed that the star's velocity is parallel to its 
rotation axis uj« . The standoff distance of the shock z s h 
is given by equation (3). The pulsar wind extends out 
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to a radial distance r m o ~ z s h from the z— axis at z = 0. 
The magnetic field strength at this distance is 



or 



B m0 = B L ( — ) = (4:Trp lsm v 2 ) 

V r rnO ) 



1/2 



4.6 x 10 



-4 1/2 
U ism 



v 8 G 



(25) 



which is independent of the period and magnetic moment 
of the pulsar. The magnetic field pressure of the wind at 
r m o is 

B 2 

r p lsm M 2 , (26) 



PmO = 



4ir 



where pi sm is the ambient pressure of the interstellar 
medium, T is the ratio of specific heats, and M. = v/c s 
is the Mach number of the pulsar. Typically, M > 1, 
so that electromagnetic pressure at r w is much larger 
than pi sm . The pressure of the shocked interstellar gas is 
larger than pi sm by a factor M 2 neglecting radiative cool- 
ing. Because the time scale t rn = r m o/v is much shorter 
than the spin-down time scale t, the time-dependence of 
the pulsar's period can be neglected. 

The momentum of the relativistic pulsar wind is de- 
flected by the internal shock as indicated in Figure 3. 
The flow of power and magnetic flux is however essen- 
tially unchanged. We assume that the internal shock 
has the effect of establishing an approximate equipar- 
tition between the field energy density and the energy 
density of relativistic particles. It is unlikely that the 
pulsar's magnetized wind is compressed into a thin layer 
in contact with the interstellar medium as discussed in 
the hydrodynamic model (Baranov, Krasnobaev, & Ku- 
likovskii 1971; Wilkin 1996). Plasma in the magnetotail 
moves away from the pulsar with an initially relativistic 
velocity, v m « v mz z, where v mz (z = 0) = 0(c). As indi- 
cated in Figure 3, the magnetic field in the magnetotail 
is predominantly toroidal, B w cj)B m (z). 

Inside the magnetotail there is a cylindrical neutral 
layer across which the magnetic field reverses direction 
(see also simulations by Romanova et al. 2001; Toropina 
et al. 2001). Reconnection or annihilation of the mag- 
netic field in this layer continually accelerates the charged 
particles of the flow. In the following equations, we 
assume that the charged particles consist of electrons 
and positrons. The flux of particles from the pulsar is 
dNi/dt w 2£>£r£c/e. In steady state this particle flux 
flows down the magnetotail so that dNg/dt = Trr m ne v mz , 
where ne is the lepton density. The equipartition Lorentz 
factor of the leptons is 

1 < Hi. ' I. ' 
< It >= Tr- 



lO 6 M3£ 



By, 



8 mc 2 \ B m0 r m0 



BrnTn 



(27) 



P 2 \B m0 r m0 J V!0 10 cm/s 

where we have used the fact that B m or m o = B^rr.. 

The total, kinetic plus field, energy flux along the mag- 
netotail is 



Fe(z) = -(B m r m ) 2 v„ 



(28) 



Assuming a stationary flow with respect to the pulsar, 
the synchrotron losses along the magnetotail imply 

dj~ E , n x 2 o t-,2 2 

-^T = -^r m n t )-r e cB m < J e > , 



d{x 2 v mz ) 1 / 

dz' ~ ~~ \~ 1 X Vmz ■ 



(29) 



where z' = — z > 0, and 

X = (B m /B m0 )(r m /r m0 ) < 1 . 
Here, zo is the radiation length-scale of the magnetotail, 

24m 2 c 4 



zo = 



B L r L r 2 e e (4Trp ism v 2 ) C e 

p2 



15 



M30 nism V g U£ 



2 p PC , 



(30) 



and Ct =< 7 2 > / < > 2 is a constant of order unity. 
Solutions of equation (29) for v mz — const, indicate that 
X decreases very gradually with z' . 

The frequency of the peak of the synchrotron radiation 
from the magnetotail is 



y syn 



1 eBml 2 
47r mc 



(31) 

where we have used equation (27) for lepton Lorentz fac- 
tor. 

Most pulsar velocities through the ISM are less than 
10 3 km/s, with typical velocities ~ 100 — 500 km/s (Ar- 
zoumanian, Chernoff, & Cordes 2002). For these slower 
pulsars, the length zq may be larger than the total dis- 
tance the pulsar has travelled. For such conditions only 
a small fraction of the magnetotail energy flux Te is lost 
to synchrotron radiation. Here we discuss the long dis- 
tance z' » r m0 behavior of the magnetotail in this limit. 
The internal pressure of the magnetotail, p m , is much 
larger than the pressure of the ISM at the distance z' , 
Pism — PismC 2 /T <C p m - In general the radius of the 
magnetotail increases with distance from the star. As- 
suming a stationary configuration in the reference frame 
moving with the star we have 



dr m (z') _ v, 



Thus 



Pm(z') 

Pis 



dr m (z') 
dz' 



dz' 



mr 



1/2 



Zltr m VmzPi. 



(32) 



4c 



1/2 



r m. Vm.: 



0, , I ''mO 
.02 v% [ 

' m. 



10 10 cm/s 



3/2 



(33) 



Expressed in this form, dr m /dz' is independent of E w 
and pism- Equation (32) represents a balance of the in- 
ternal pressure of the magnetotail, p m (z'), against the 
ram pressure due to the tail's expansion into the ISM, 
p lS m(dr m /dt) 2 , where dz' /dt = v mz - 
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Fig. 4. — Qualitative dependence of the width of the magnetotail r m (z') on the distance behind the star z'. The vertical scale is arbitrary; 
it given by (z' / L) 1 / 2 F(z' / L), with F given below equation (34). 



For distances z' ^> r m Q but small compared with the 
total distance the star has travelled, v mz is approxi- 
mately constant. Consequently, 

. m (z0«O.2(,„ lO .0 1/2 f^^) 3/4 , 04) 

\ "mz / 

for z' > r m0 . 

The situation is different for distances z' of the order of 
or somewhat larger than the total distance travelled by 
the star, L = vt « 10ws(</10 4 yr) pc, with t the age of the 
star. The velocity of the magnetotail plasma v mz must 
decrease strongly with z' of the order of this distance ow- 
ing to collisions with the plasma of the above-mentioned 
bubble TZ(t) formed at the birth of the pulsar. For a 
rough description consider 



v m (z') 



»(0)exp - 



A2 



(il 

L 2 



For this case equation (33) gives 
r m (z) = 0.2(r m0 ,) (^J^) F [j ) , 



where 



(35) 



F®= (j^exp(3£ 2 /2) 



1/2 



Figure 4 shows the qualitative dependence of r m (z'). 
Notice that the curvature of the magnetotail radius 
d 2 r m /dz' 2 is initially negative but at large z' > 0.45L 
it becomes positive. Further note that v mz (z / )[r m (z')] 2 
is an increasing function so that B m (z') decreases with 
distance from the star. 

The magnetotails may be observable as low sur- 
face brightness regions of non-thermal, polarized (syn- 
chrotron) emission. 



6. EXAMPLES OF MAGNETOTAILS: GUITAR AND MOUSE 
NEBULAE 

The Guitar Nebula is created by the high velocity 
pulsar B2224+65. Chatterjee and Cordes (2002, 2004) 
estimate the velocity of the pulsar as vg « 1.7 or 
1700 km/s based on a distance to the source of D = 
1.9 kpc. The pulsar period P — 0.68 s and spin- 
down rate P = 9.7 x 10~ 15 s/s imply a power output 
E w 1.2 x 10 33 erg/s assuming a moment of inertia 
I = 10 45 g cm 2 . The angular length of the tail is 
69 w 15 (Chatterjee & Cordes 2004) which corresponds 
to a tail length z 1 w 0.15 pc/sint, where i is the an- 
gle between the jet axis and the line of sight. With 
E = B 2 L r 2 L c, this power implies ^30 ~ 2.1. Consequently, 
equation (30) gives the length-scale of the magnetotail 
z ~ 1 pc/(n ism Ce). This corresponds approximately to 
the whole length of the tail in the Guitar nebula ob- 
served in the H a line. Equation (31) gives the estimate 

v syn « 2.6 x 10 16 Hz n\i 2 m x 5 (r m0 /r m )(v m z/10 10 cm/s) 2 . 
The Guitar nebula was observed as a source of very weak 
radiation in soft X-rays which is a sign of the presence 
of highly relativistic leptons, 7 ~ 10 7 . (Romani, Cordes, 

6 Yadigaroglu 1997). Equation (13) also gives similar 

7 for the Guitar nebula. Both, the magnetic field and 
the particles interact with the interstellar medium and 
are deflected to the tail. The radio images of the Gui- 
tar nebula show clear limb- brightening (e.g., Chatterjee 
& Cordes 2004). This may be due to the compression 
of the magnetic field at the shock wave which leads to 
higher intensity of radiation (e.g., Romani et al. 1997). 
From other side, enhanced reconnection is expected near 
the shock wave due to compression the magnetic field 
components of opposite polarity. This will also lead to 
limb-brightening. The inner (smaller scale) structure of 
the Guitar nebula is closed in the back and may repre- 
sent the inner shock wave, where accelerated particles 
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were stopped by the interstellar medium. If this is the 
case, then the flow has a < 1 before the shock wave. 

Another example is a Mouse Nebula which is created 
by much more powerful pulsar PSR J1747-2958 which 
propagates with high velocity v « 600 km/s (Camilo et 
al. 2002). It has a period P = 0.098 s and a spin- 
down luminosity E k, 2.5 x 10 36 erg/s which corresponds 
to ^30 ~ 2.5. The extended radio nebulae was discov- 
ered around this pulsar (Yusef-Zadeh & Bally 1987). 
An interesting feature is the tail which propagates to 
very far distances. Recently an X-ray nebula was found 
in observations with Chandra (Gaensler et al. 2004). 
Gaensler et al. (2004) estimate the X-ray luminosity of 
the nebula as L x (0.5 - 8 keV) « 5 x 10 34 erg/s based 
on a distance to the source of D — 5 kpc, the angu- 
lar length of the Mouse tail as 59 « 45 which corre- 
sponds to a tail length z' ~ 1 pc/sint. For these pa- 
rameters and assuming v$ = 0.6 and rii sm = 0.3/cm 3 
(Gaensler et al. 2004), equation (30) gives radiation 
length-scale of the magnetotail, zq w 0.5 pc/C^. For 
this object equation (31) gives the estimate iy syn ~ 2.7 x 

10 19 Hz nV s l x 5 (Wr m )(tWl0 lo cm/s) 2 . The pre- 
dieted frequency and overall luminosity is much higher 
than those for Guitar nebula, which correspond to ob- 
servations of much brighter X-ray source. 

The shape of the Mouse nebula was modelled by hy- 
drodynamic simulations (Gaensler et al. 2004; van der 
Swaluw et al. 2003; and Bucciantini 2002). It was sug- 
gested that the inner brightest region of the nebula may 
be due to the termination shock of the pulsar wind, while 
the long tail may represent the post-shock flow. We 
note that the Chandra observations do now show evi- 
dence of the shock wave in the X-ray picture of this neb- 
ula. Instead, the X-ray luminosity in the brightest region 
(the head of the Mouse) varies gradually (Gaensler et al. 
2004) . This may be a sign that in the tail the a parame- 
ters is not small and the back flow is not super- Alfvenic. 
If a ~ 1 in the tail, then the magnetotail may propagate 
to very far distances (see, e.g., simulations by Romanova 
et al. (2001) and Toropina et al. (2001). Radio radi- 
ation is then connected with particles which propagate 
from the pulsar along the tail to far distances and lose 
their energy there, or, particles may be accelerated in 
situ during the reconnection processes in the magneto- 
tail (Toropina et al. 2001). 

The direction of jets observed in a few cases is almost 
aligned with the direction of propagation of the pulsar 
and seems to be perpendicular to the PWN disk (e.g., 
Ng & Romani 1987). If the direction of the pulsar's 
motion is determined by a magnetic kick during its for- 
mation (Ardeljan, Bisnovatyi-Kogan & Moiseenko 2001; 
Lai, Chernoff & Cordes 2000) then both, u>* and fi will 
be in the same direction. If the magnetic axis is mis- 
aligned relative to the rotational axis, but w» is aligned 
with the pulsar's velocity v, then the situation is similar 
to the one described above, because the jet forms in the 
direction of w* and the disk magnetic structure is per- 
pendicular to w*. However, if u>* does not coincide with 
v, then both the jet and the PWN disk may give non- 
axisymmetric features of the pulsar wind nebula. This 
later situation has not been investigated. A related sit- 
uation appears in the recently discovered double pulsar 
binary PSR J0737-3039 where the direction of the flow 



from the wind-generating pulsar changes in time because 
of the motion of stars in the binary system (Kaspi et al. 
2004; Demorest et al. 2004). Three dimensional simula- 
tions of such a flow have shown that for large inclination 
angles between magnetic and rotational axes of pulsar, 
the shock wave may be non-axisymmetric (Spitkovsky & 
Arons 2004). 

7. OBSERVABILITY OF WIND AND MAGNETOTAIL 
EMITTING STARS IN THE SOLAR NEIGHBORHOOD 

Although neutron stars which spin-down to periods 
P > 1 — 3 s shut off as pulsars, they may continue to 
emit a wind and interact with the ISM. The number of 
such objects is larger than the number of short period 
pulsars. It is possible that such objects are detectable at 
distances < 10 2 pc. Interesting nearby isolated neutron 
star candidate RX J1856. 5-3754 has a prominent shock 
wave, and we discuss this object below. 

Different spin-down powers E w are expected depending 
on the star's period of rotation and magnetic field. We 
consider two cases, one with a relatively high spin-down 
power E w ~ 10 31 — 10 32 ergs/s which is pertinent to old 
pulsars, or relatively powerful shutoff pulsars, and a low 
spin-down power, E w — 10 28 ergs/s corresponding to very 
weak shutoff pulsars. 

7.1. Relatively High Spin-Down Power Objects 

Shutoff or misdirected pulsars with periods P ~ 1 s 
and a surface magnetic field B ~ 10 12 G have a rela- 
tively high spin-down power, E w ~ 6 x 10 31 erg/s (i.e., 
equation 2). The synchrotron emission from the pulsar 
wind for distances within the bow shock wave is given by 
equation (20) and the photon energies are given by equa- 
tion (22) and are in the high X-ray range. Additionally, 
there will be a small power in thermal radiation from the 
shocked ISM according to equation (5). For high veloc- 
ity pulsars, say v > 500 km/s, a significant fraction of 
the spin-down power is emitted in the magnetotail. The 
photon energies are then much lower, in the UV band 
according to equation (31). Such shut-off or misdirected 
pulsars may be observed owing to the radiation of their 
tails. 

7.2. Nature of RX J 1856.5-3754 

The source RX J1856. 5-3754 is an interesting isolated 
neutron star candidate located at a distance d « 60 pc 
(Walter 2001). It has a bow shock wave observed by 
van Kerkwijk and Kulkarni (2001a, b) in the H Q line. 
The observed standoff distance between the star and the 
bow shock wave is « 1" which corresponds to z s h « 
0.9 x 10 15 cm. The bow shock is probably due to the in- 
teraction of star's wind with the ISM (see, e.g., van Kerk- 
wijk and Kulkarni 2001b). In this case we obtain from 

equation (3) the relation Bi 2 / {n^^v^P 2 ) « 9, where B\ 2 
may be replaced by /130. This relation implies that the 
spin-down power of the star is E w w 5.2 x 10 32 n ism v 2 
erg/s w 5.2 x 10 30 erg/s, where we took w 8 = 0.1 and 
n ism = 1/cm 3 . 

Clearly, a range of values of B 12 (or ^30), P , and 
v 8 can give ^3o/(^8-P 2 ) ~ 9, where we have taken for 
simplicity ni sm = 1 (1/cm 3 ). To consider the likely 
values, note that the spin-down time of the star is 
t sd = P/P = /P 2 c 3 /(47rV) « 2.2 x 10V P 2 /l4o, 
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assuming I = 10 g cm 2 . The total number of neu- 
tron stars likely to be within a distance of 60 pc of the 
Sun is N w 10' 9 (60pc/10 4 pc) 3 « 220 (e.g., Schvartsman 
1978; Treves & Colpi 1991). Then the probable num- 
ber of objects with spin-down time t sd within 60 pc is 
N' w 220(t ad /t H ), where t H ~ 1.5 x 10 10 years is the 
Hubble time. Thus, for the object RX J1856.5-3754, 
N' w 0.036/(/X3of8)- The corresponding rotation period 
of the star is P w 0.33(^3o/ws) 1 ^ 2 s. Clearly, larger values 
of AT' occur for a smaller values of the product ^.30 fs • For 
example, for /130 = 0.1 and vs = 0.1, N' « 3.6. For these 
values E w w 5.2 x 10 30 erg/s. So, this object may be a 
misdirected pulsar if magnetic field is large, or, shutoff 
pulsar if the field is small. 

7.3. Low Spin-Down Power Objects 

There should be a much larger number of shutoff pul- 
sars with the lower spin-down power. As an example, 
we take objects with the spin-down power E w = 10 28 

erg/s. For these objects, /z| /P 4 w 1.7 x 10~ 4 , so that 
1/2 

P « 8.7/i 3 Q s. The number of such objects within a dis- 
tance of 60 pc of the Sun is AT' « 220(t sd /t H ) w 24//z 30 . 
The standoff distance of the bow shock wave is however 
small, z s h w 1.3 x 10 12 Wg -1 cm, but much larger than the 
light cylinder radius. 

8. DISCUSSIONS AND CONCLUSIONS 

This work considers rotating magnetized neutron stars 
emitting power E w in a wind of relativistic particles and 
electromagnetic field. We argue that a non-negligible 
fraction of the wind's power may be converted to rela- 
tivistic particles due to annihilation or reconnection of 
the magnetic field. Outside of the light cylinder, the 
star's rotation acts to wind up the magnetic field to form 
a predominantly azimuthal, slowly decreasing with dis- 
tance, magnetic field of opposite polarity on cither side 
of the equatorial plane normal to the star's rotation axis. 
An analogous situation exists in the Solar Wind (Bcrtin 
& Coppi 1985). The magnetic field annihilates across the 
equatorial plane with the magnetic energy going to ac- 
celerate the charged particles to highly relativistic ener- 
gies. The accelerated leptons emit synchrotron radiation 
in a broad range from the UV to gamma ray energies. 
Additionally, there may be oppositely directed Poynting 
outflows along the star's rotation axis. These outflows 
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are analogous to the Poynting jets discussed by Lovelace 
and Romanova (2003). The model is qualitative in the 
respect that we do not have global self-consistent calcu- 
lations of the electromagnetic fields and sources. 

For a typical, supersonically moving star, the star's 
relativistic wind forms a bow shock wave with the in- 
terstellar medium (Schvartsman 1970). We argue that 

an appreciable fraction of the star's wind power E w is 
deflected by the bow shock wave and collimated into 
the star's magnetotail. Plasma moves down the mag- 
netotail with a relativistic velocity. Further, we argue 
that equipartition is reached in the magnetotail between 
the magnetic energy and the relativistic particle energy. 
For the case where the charged particles are leptons, the 
synchrotron radiation length-scale zq of the magnetotail 
is calculated. Over this distance the energy flux in the 
magnetotail decreases by a factor of order 2. For highly 
supersonic pulsars, zq may be less than the total distance 
the star has travelled. The synchrotron radiation spec- 
trum is expected to be nonthermal with typical photon 
energies in the UV range for the highly supersonic pul- 
sar B2224+65 which generates a bow shock wave and the 
Guitar Nebula (Chatterjee & Cordes 2002). 

The ratio of the power due to the shock heating of the 
ISM, which gives thermal emission, to the wind power 
E w , which gives nonthermal synchrotron emission, is 
shown to be (3/64) (v/c) <C 1, where v is the velocity 
of the star and c the speed of light. 

An equation is derived for the radius of the magnetotail 
r m {z') as a function of distance z' from the star. For large 
distances z', of the order of the distance travelled by the 
star, we argue that the magnetotail has a 'trumpet' shape 
owing to the slowing down of the magnetotail flow. 

We argue that the isolated neutron star candidate RX 
J1956. 5-3754 may be a misdirected or shutoff pulsar. 

We estimate the number of shutoff pulsars which may 
be observable in the vicinity of the Sun for cases of rela- 
tively strong (10 32 erg/s) and weak (10 28 erg/s) power. A 
much larger number of weak shutoff pulsars is predicted. 
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